Energy Gaps and Stark Effect in Boron Nitride Nanoribbons 
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A first-principles investigation of the electronic properties of boron nitride nanoribbons (BNNRs) 
having either armchair or zigzag shaped edges passivated by hydrogen with widths up to 10 nm 
is presented. Band gaps of armchair BNNRs exhibit family-dependent oscillations as the width 
increases and, for ribbons wider than 3 nm, converge to a constant value that is 0.02 eV smaller 
than the bulk band gap of a boron nitride sheet owing to the existence of very weak edge states. The 
band gap of zigzag BNNRs monotonically decreases and converges to a gap that is 0.7 eV smaller 
than the bulk gap due to the presence of strong edge states. When a transverse electric field is 
applied, the band gaps of armchair BNNRs decrease monotonically with the field strength. For the 
zigzag BNNRs, however, the band gaps and the carrier effective masses either increase or decrease 
depending on the direction and the strength of the field. 



Two-dimensional crystals, including graphene and sin- 
gle layer of hexagonal boron nitride (BN), have recently 
been fabricated [l|. Among them, only graphene has 
been studied extensively [2]. Unlike graphene, a hexago- 
nal BN sheet is a wide gap insulator like bulk hexagonal 
BN jsf and is a promising material in optics and opto- 
electronics y. 

Graphene nanoribbons (GNRs) [5] with width a few 
to a hundred nanometers have been produced by litho- 
graphical patterning 0] or chemical processing [8] of 
graphene. We expect that boron nitride nanoribbons 
(BNNRs) could also be made using similar or other tech- 
niques. Figures [Tfa) and [Tfb) show the structures of 
an armchair BNNR with Ng, dimer lines (A/'a-aBNNR) 
and a zigzag BNNR with zigzag chains (A/'z-zBNNR), 
respectively. A tight-binding study of the bandstruc- 
tures of 21-aBNNR and 13-zBNNR (corresponding to 
widths ~ 3 nm) [9] and first-principles investigations of 
the electronic properties of small width BNNRs 0, |Tl[ 
have been reported. However, to our knowledge, first- 
principles calculations on the electronic properties of ex- 
perimentally realizable size of BNNRs have not been per- 
formed. 

Under a transverse electric field, carbon nanotubes 
with impurity atoms are expected to show novel band gap 
opening behaviors [12], whereas zigzag GNRs reveal half- 
metallicity On the other hand, single- walled boron 
nitride nanotubes (SW-BNNTs), which are rolled up BN 
sheets [T3, [lH, [l6| , have been predicted to show gigantic 
Stark effect in their band gaps in response to a transverse 
electric field [171, and this effect has been confirmed ex- 
perimentally [18]. The effect becomes stronger in larger 
diameter SW-BNNTs [l^- A similar phenomenon is ex- 
pected in BNNRs. Unlike SW-BNNTs, however, BNNRs 
can be arbitrarily wide. Therefore, the consequences of 
the Stark effect in BNNRs would be even more dramatic 
than in SW-BNNTs. 

In this study, we report first-principles calculations on 
the electronic properties of armchair and zigzag BNNRs 
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FIG. 1: Schematic of (a) 14-aBNNR and (b) 7-zBNNR pas- 
sivated by hydrogen atoms. Boron, nitrogen and hydrogen 
atoms are represented by white, black and grey spheres, re- 
spectively. BNNRs are periodic along the y direction. 



up to width of 10 nm with hydrogen passivation of the 
edge carbon atoms. We show that the band gaps of arm- 
chair and zigzag BNNRs do not converge to the same 
value even when the ribbons are very wide. The band 
gap of armchair BNNRs, obtained by density functional 
theory (DFT) calculations within the local density ap- 
proximation (LDA), converges to a value that is 0.02 eV 
smaller than the LDA band gap of 4.53 eV of a BN 
sheet [19]. Unlike armchair GNRs, the lowest unoccu- 
pied band of the armchair BNNRs is composed of edge 
states with energy position asymptote to a fixed value 
when the ribbon is wider than 3 nm, the decay length of 
the edge-state. The band gap of the zigzag BNNRs, also 
determined by edge states, is monotonically reduced as a 
function of increasing width and converges to a value that 
is 0.7 eV smaller than the LDA bulk gap because, as dis- 
cussed below, of an additional edge polarization charge 
effect. The DFT Kohn-Sham eigenvalues within LDA in 
general underestimate the band gaps of materials; an ac- 
curate first-principles calculation of band gaps requires a 
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quasiparticle approach [20| . The basic physics discovered 
here however should not be changed. 

When a transverse electric field is applied, the highest 
occupied and the lowest empty states in armchair BNNRs 
become localized at the two different edges. Because of 
the external electrostatic potential difference between the 
two edges, the band gap is reduced with increased field 
strength. On the contrary, in zigzag BNNRs, depending 
on the field direction, the states near the band gap either 
become more localized at the edges or less so. Also, the 
band gaps and effective masses either decrease or increase 
depending on field strength and direction. These novel 
properties could be used in manipulating the transport 
properties of doped BNNRs. 

We performed ab initio pseudopotential DFT calcula- 
tions within LDA in a supercell configuration using the 
SIESTA computer code [21,] . A double-zeta plus polariza- 
tion basis set was used and ghost orbitals were included 
to describe free-electron-like states BE 111. A charge 
density cutoff of 300 Ry was used and atomic positions 
were relaxed so that the force on each atom is less than 
0.04 eV/A. To eliminate spurious interactions between 
periodic images, a supercell size of up to 20 nmx20 nm 
in the xz plane was used. 

The armchair BNNRs are found to have a direct gap 
at the zone center [left panel of Fig. [2l^a)]. The high- 
est occupied state, the valence band maximum (VBM), 
has a wavefunction which is localized at nitrogen atoms 
throughout the ribbon [right lower panel of Fig. [2l^a)]. 
The lowest empty state, the conduction band minimum 
(CBM), is however an edge state with wavefunction local- 
izing at the boron atoms on the edges [right upper panel 
of Fig. [21(a)]. In contrast, the corresponding state for an 
armchair GNRs is delocalized throughout the ribbon [23]. 
The total potential near the edge of the armchair BNNRs 
is different from that of the bulk By incorporating 
this variation into the on-site potential energies of a few 
BN dimer lines near the edges, one could reproduce the 
main features of the states near the band gap within a 
tight-binding formulation [24]. 

The zigzag BNNRs have the VBM at a point between 
the X and the F points which has wavefunction localized 
at the nitrogen edge and the CBM at the X point which 
has wavefunction localized at the boron edge [Fig. [31(a)], 
in agreement with results by Nakamura et al [10] 

Figure [Ha) shows the band gap variation of the arm- 
chair and the zigzag BNNRs with width. The most in- 
teresting and somewhat counter-intuitive feature is that 
as the width increases, band gaps of the armchair and 
the zigzag BNNRs converge to different values neither of 
which is that of a BN sheet. This is because the CBM or 
the VBM are determined by edge- states. 

The edge-state band gaps of armchair BNNRs show a 
family behavior with respect to the number of dimer lines 
A^a [inset of Fig.[3Ka)], the family with A^a = 3n — 1 having 
the smallest gaps where n is an integer. A similar trend 
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FIG. 2: (color online) LDA energy bandstructure (left) and 
the squared wavefunctions integrated along z of the highest 
occupied state (right lower) and the lowest unoccupied state 
(right upper) of 14-aBNNR under an external electric field 
^ext of strength (a) zero and (b) 0.1 eV/A along +x direction. 
Dashed red lines in the bandstructure indicate the energies of 
the band edge states. In the wavefunction plots, green regions 
are associated with high densities. 
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is observed in armchair GNRs 
edge-state band gaps of the armchair BNNRs converge to 
a near constant value roughly when the ribbon is wider 
than 3 nm [Fig.[ll^a)]. This characteristic length is related 
to the decay length of the CBM edge-state. Figure [H^b) 
shows the squared electron wavefunctions of the CBM 
states of 14-aBNNR and 26-aBNNR integrated in the yz 
plane. These states are localized on the boron atoms near 
the two edges. When the width is about 3 nm as in the 
26-aBNNR, the wavefunction from the two edges begins 
to decouple and thus stabling its energy position. 

In zigzag BNNRs, the boron edge and the nitrogen 
edge are negatively and positively charged (electronic 
plus ionic charge), respectively. Because of this polar- 
ization, the potential felt by electrons is higher at the 
boron edge and lower at the nitrogen edge, contribut- 
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FIG. 3: (color online) LDA energy bandstructure (left) and 
the squared wavefunctions integrated along z of the highest 
occupied state (right lower) and the lowest unoccupied state 
(right upper) of 7-zBNNR under an external electric field ^ext 
of strength (a) zero, (b) 0.1 eV/A and (c) -0.1 eV/A along the 
X direction. Dashed red lines in the bandstructure indicate 
the energies of the band edge states. In the wavefunction 
plots, green regions are associated with high densities. 



ing a factor which increases the band gap of the narrow 
zBNNRs since the VBM and the CBM are edge-states 
localized at the nitrogen edge and the boron edge, re- 
spectively [Fig. [3l^a)]. However, as the ribbon becomes 
wider, the effective polarization line charge density (Jeff, 
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FIG. 4: (color online) (a) Band gaps of armchair (filled red 
squares) and zigzag (empty blue squares) BNNRs versus their 
widths. Inset: band gaps of armchair BNNRs versus A^a (see 
Fig. [1]). Solid lines are a guide to the eyes. Dashed lines 
indicate the bulk band gap of a BN sheet with no edges, (b) 
Probability distributions |$(r)|^ integrated in the yz plane 
(see Fig. [1]) versus the distance along the x direction from the 
ribbon center for the lowest unoccupied state in 14-aBNNR 
(soUd red hne) and 26-aBNNR (dashed blue line), (c) The 
effective polarized line charge density (Jeff of zigzag BNNRs 
versus Wz- The solid line is a guide to the eyes. 
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where p{r) is the total charge density including the core 
charge, and the spatial period along the y direction 
[see Fig. [TJb)], decreases as ~ l/it;^ [Fig. [H^c)] due to an 
increased screening, resulting in the decrease and conver- 
gence of the band gap as Wz increases. 

Figure [21(b) shows how the bandstructure and wave- 
functions of a 14-aBNNR change under a 0.1 eV/A trans- 
verse electric field. Owing to the Stark effect, the wave- 
functions of the highest occupied and the lowest unoc- 
cupied states now localize at the opposite edges where 
the external electrostatic potential felt by an electron is 
higher and lower, respectively [right panel of Fig. EJb)]. 
Thus, the band gaps of armchair BNNRs decrease when a 
transverse electric field is applied [left panel of Fig. [21(b)]. 
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FIG. 5: (color online) (a) and (c): LDA energy gaps and effec- 
tive hole masses (in units of the free electron mass mo) of the 
highest occupied band in 36-aBNNR (filled red circles) and 
84-aBNNR (empty blue squares) under a transverse electric 
field versus the field strength. The inset in (a) shows energy 
gaps as a function of the external potential difference between 
the two edges, (b) and (d): Similar quantities as in (a) and 
(c) for 27-zBNNR (filled red circles) and 46-zBNNR (empty 
blue squares), respectively. 



A similar behavior has been predicted [IT] and ob- 
served [18] in BNNTs. 

Figures [3](b) and[3](c) show how the bandstructure and 
the edge-state wavefunctions of a 7-zBNNR change un- 
der a 0.1 eV/A transverse electric field. When an elec- 
tric field is applied toward -\-x direction, the VBM and 
CBM edge-state wavefunctions do not change qualita- 
tively [right panel of Fig. [Sjb)]. Thus, for a similar rea- 
son as in the armchair BNNRs, the band gap decreases 
[left panel of Fig. [3](b)]. When the electric field is ap- 
plied along —X direction, the potential felt by an elec- 
tron localized at the right edge (the nitrogen edge) is 
decreased whereas that at the left edge (the boron edge) 
is increased. Therefore, the energy gap between these 
two states increases as shown in Fig. [3fc). (Actually, 
the energy of the original lowest empty state has been 
pushed upward by so much that the state is no longer 
the CBM state.) At the same time, the VBM states tend 
to delocalize [right lower panel of Fig. ^c)]. 

Figures [Sl^a) and [S^b) show the band gap variation of 
BNNRs with field strength. In armchair BNNRs, the 
band gap decreases when the field strength increases re- 
gardless of its direction [Fig. El^a)]. A similar behavior 
has been observed in SW-BNNTs [17]. For example, for 
the 10 nm wide 84-aBNNR the LDA band gap is reduced 
from 4.5 eV to less than 1.0 eV under a 0.1 eV/A field. In 



zigzag BNNRs, at small field strength, the band gap de- 
creases when the field is along -\-x direction but increases 
when the field is reversed [Fig. [Sfb)]. In other words, 
zigzag BNNRs show asymmetric Stark effect. However, 
as the field becomes stronger, the band gap decreases re- 
gardless of the direction. Moreover, the band gap varia- 
tions, when plotted against the difference in the external 
electrostatic potential between the two edges, fall on a 
universal curve for ribbons with different widths [insets 
of Figs. [5](a) and[5](b)]. This is because the gap deter- 
mining states localize on different edges as the field be- 
comes strong; thus, the change in their energy difference 
is directly related to the potential difference between the 
edges. 

Figures [5jc) and ^d) show the effective mass of the 
hole carrier at the VBM for a range of external field 
strength. In armchair BNNRs, the hole mass of the VBM 
is independent of the field strength. In contrast, the 
corresponding effective mass in zigzag BNNRs changes 
with the external field, and even more interestingly, in 
an asymmetric way. In particular, when the field is 
along —X direction, the effective mass decreases substan- 
tially. Within ±0.02 eV/A variation of the field, the 
effective mass can be varied by 50 % from 0.6 mo to 
0.9 mo where mo is the free electron mass. In the case of 
electron carriers in zigzag BNNRs, a nearly-free-electron 
state j3, US, 0] becomes the CBM if a field stronger 
than a critical value, depending on the width, is applied 
[see Fig. Efc)], and the characteristics of charge carri- 
ers change significantly. These novel phenomena demon- 
strate the possibility of tuning carrier mobilities of doped 
BNNRs by applying a transverse electric field. 

In summary, we have studied the electronic proper- 
ties of BNNRs as a function of width with or without 
an external transverse electric fields. The band gap of 
the armchair BNNR and that of the zigzag BNNR are 
determined by edge-states and thus converge to values 
different from that of the bulk BN sheet. The electronic 
and the transport properties of BNNRs are shown to be 
tunable by an external transverse electric field. Espe- 
cially, zigzag BNNRs are shown to exhibit asymmetric 
response to the electric field. 

Additional remark: After completion of the work and 
during the preparation of this manuscript, we became 
aware of a related work on similar systems from other 
group [29]. 
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